OHSTPY-HEP-T-96-029 
September 1996 



PRODUCTION OF THE B, MESON* 



YU-QI CHEN 

Department of Physics, Ohio State University 
Columbus, OH 43210, USA 

The production of the Be and B* mesons was studied in the framework of the factor- 
ization formahsm and perturbative QCD. Predicted results imphed that an observable 
number of Be and B* events can be produced at LEP I and at Tevatron. The fragmen- 
tation approximation describes the production of the B^ and B* mesons very well in 
high energy e+e" collisions. In hadronic collisions, it is valid when and only when the 
transverse momentum Py of the produced Be and B* is much larger than the mass of 
the Be meson. 

1. Introduction 

The Be meson is very interesting because it is the only flavored meson containing 
a heavy quark and a heavy antiquark. However, the study of the Be meson suffers 
from the difficulty of producing it in experiments. The conventional production 
mechanisms for the Bu and Bd mesons and for heavy quarkonia J/ip and T are 
not available for the Be and B* mesons. A question is whether a sufficiently large 
number of Be events can be produced in experiments. If not, the study of the Be 
meson would be only of theoretical interest. 

Theoretical studies of the production for the Be meson over the past few years 
give a positive answer to this question. The dominant mechanism for the production 
of the Be meson in high energy collisions has been found to be the following. First 
two heavy quark pairs cc and bb are produced in a high energy process, and then the 
c quark from the cc pair and the b quark from the bb pair bind to form a Be mesc!iL| 
At the e+e^ collider LEP 1, the Be meson can be produced via the decayu'n 
e+e~ — > Z'^ ^ i3c + 6 + c. At the hadron-hadron colliders Tevatron and LHC, it can 
be produced by gluon-gluon fusion and quark-antiquark annihilation subprocessesU 
such as g + g Be + b + c, q + q ^ Be + b + c. Calculations show that quite a large 
number of the Be events can be produced both at LEP I and at the Tevatron and 
the LHC. The experimental search for the Be meson is now under way at LEP and 
the Tevatron. Some preliminary results from LEP I experiments were reported in 
this workshops. Results from the Tevatron are expected in the near future. 

Like the cc and bb systems, the be is a nonrelativistic bound state system of 
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a heavy quark and a heavy antiauark. According to the nonrelativistic QCD 
(NRQCD) factorization formahsmO, the production cross section of the Be and 
the B* mesons can be factorized as the products of short distance coefficients and 
long distance matrix elements. The production of the cc and the bb pairs is a short- 
distance process that occurs at a distance scale of order 1/toq or smaller, and it 
can be accurately calculated using perturbative QCD to any order in as(mQ). The 
formation of the Be meson is a long-distance process which occurs at a scale of order 
of l/{mQv), and it can be described by NRQCD matrix elements which scale in a 
definite way with v, the typical relative velocity of the bound states. In the produc- 
tion of the Be or B* mesons, the contribution from color-octet processes is smaller 
than that from color-singlet process simply because the short-distance coefficients 
are of the same order in as , but the color-octet matrix elements are higher order of 
V. Thus the production is dominated by color-singlet processes. The relevant color- 
singlet matrix elements can be estimated well by heavy-quark potential models or 
by lattice calculations. Therefore, the production of the Be is predictable without 
any free parameters. 

If the Be is produced with sufficiently large transverse momentum, then the 
fragmentation approximation is expected to be available. In this approximation, 
the cross section factors into a cross section for producting a b quark and a frag- 
mentation function that gives the probability for the b quark to hadronize into a Be 
meson. The fragmentation functions arc universal, so they can be used to calculate 
the production of the Be both in e+e^ process and in hadron-hadron collisions. 
An interesting question is how well the fragmentation approximation to the full 
calculation works for the production of the Be and B* mesons. In the following 
sections, this question will be addressed for the production of Be and B* at LEP 
and at hadronic colliders. 

2. Production of the Be Meson in e+e Collision 

At LEP I, millions of Z° events have been accumulated. Such a large number of 
events provides us with an opportunity to search for the Be meson through the 
process Z° BeiB*) + b + c. Calculations of the production rateOoB show that 
the branching ratio for Z'^ ^ Be + X is around 10~^. 

The process can be calculated either using the covariant projection methodS or 
using the recently developed threshold expansion methoda. At the lowest order in 
as , there are four Feynman diagrams responsible for the process Z'^ Be + b + c. 
Let k, P, gi, and (72 be the momenta of the Z''^ boson. Be meson, c quark and b 
quark, respectively. We can define three Lorentz invariant variables 

_ 2fc ■ gi _ 2fc ■ 92 _ 2k-P 

which are twice the energy fractions of the Be, c quark and b quark in the Z*^ rest 
frame, respectively. They are constrained hy x + y + z — 2 and < x,y, z < I. 
The full calculation for Z^ ^ Be + X is straightforward. The decay width 
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and various distributions in x, y, and z can be expressed as products of short- 
distance coefficients and long-distance color-sing let matrix elements (0|Of (^S'o)|0) 
and {OlOi" (■^S'i)|0). Using the approximate heavy-quark spin symmetry of NRQCD, 
these two matrix elements are equal to nic and 3/^ nic, respectively, up to cor- 
rections of order w^. 

The fragmentation functions for the production of the Be can be extracted by 
taking the limit Mz/ruB^ oo. In this limit, the square of the matrix element is 
singular when the produced Be is collinear either with the b quark or with the c 
quark. In this region, either x — *■ 1 or y — > 1 but z varies from to 1. In the limit, 
the z distribution comes from those most singular terms and can be factorized as 

= - bh) ■ D-^_^^,,, (z) + r(z ^ cc) ■ D^_^^^, (z), (2) 

where D^^g(,){z) and D ^^^(t){z) are the fragmentation functions for 6 or c into a 
Be meson. The b quark fragmentation functions are 

D- (z) 4a.(2mJ^(0|Of-ego)|0) rzjl-z^ ( 

-74r -I- 68r^)z^ - 2(1 - r)(6 - 19r + 18r'^)z^ + 3(1 - rfil - 2r + 2r'^)z'^^ , (3) 
for the Be meson and 

D-b^B.iz) = ^ (l-(l^r)z)6 (^2-2(3-2r)z + 3(3 

-2r + 4r^)z^ - 2(1 - r)(4 -r + 2r^)zJ + (1 - rf(3 - 2r + 2r^)z'^ ) , (4) 



for the B* meson, where r = me/ {mi, -f m,e). The c quark fragmentation functions 
can be obtained by exchanging and me- These fragmentation functions describe 
the probability for a jet initiated by a high energy h quark to include the hadron Be 
carrying a fraction z of the jet momentum. They were first calculated by Chang and 
ChenB and Braaten, Cheung, and Yuanll. These are the lowest order fragmentation 
functions at the energy scale around toBc- All leading logarithmic terms coming 
from the collinear emission of gluons can be summed up by solving the Altarelli- 
Parisi evolution equations from the energy scale at order Mz to that at order ms^ 

d n ^ ^ ^ f'dy 



dfx 



f dy 

A J z y 



where Pi-,j (x, fj,) is the Altarelli-Parisi function for the splitting of the parton of 
type i into a parton of type j with longitudinal momentum fraction x. 
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Explicit calculations show that the result predicted by the fragmentation ap- 
proximation Eq. (||) is very close to that from the full calculationlUn. 

Integrating over the variable z, in eqs. (^); (^), the probabilities for a high energy 
b antiquark or c quark fragmentating into a Be or B* meson can be obtained. Taking 
the parameters adopted by Braaten, Cheung, and Yuann, the numerical values for 
the fragmentation probabilities are 2.2 x 10~^ for b ^ Be and 3.1 x 10~^ for b ^ B*. 
The fragmentation probabilities for c are smaller than those for b by almost two 
orders of magnitude, because it is much harder to creat an extra bb pair by emitting 
a virtual gluon than to creat an extra cc pair. 



3. Production of the Be Meson in Hadron-Hadron Collision 

Hadronic production, as first pointed out by Chang and Chenl, is dominated at 
high energies by the subprocess g + g ^ Bc{B*) + b + c. It can be calculated fully to 
order in PQCD. Since a difficult numerical calculation is involved, confirmatio. 
from some other independent calculations arc necessary. Since Chang and Chenl 
first presented the numerical results for the l-i ad ron ic prodii ction ^ calculations have 



also been done by severaLother authorsMi^ iyi2|,[L3|jl4| jl5|,g^ ^ ^j? -^^rhich are in 

agreement^ Slabospitskyll3 claimed an order of magnitude larger result than CharLg 
and ChenU. Berezhnoy et a^-ui obtained a result larger than Chang and ChenH 
and smaller than SlabospitskyH^. Masetti and Sartogolla found a result similar to 

Berezhnoy et aZ0. More recently, Berezhnoy et a/.0 found that a color factor - was 

overlooked in their previous work0. After including_this factor, their revised result 
is in agreement with Chang and ChfinU. Baranovll^ independently also obtained 
results similar to Chang and ChenB. Kolodziej et a/Jl3 presented results using 
different narton distribution functions and energy scale from that used by Chang 
and ChenU so it difficult to directly compare their final results. However their 
results for the cross sections for the parton subprocess are similar to othersEjiij. 
Kolodziej et alS^ nmnted out that the cross section for the subprocess calculated 
by Berezhnoy et alSS for y/l — 1 TeV was incorrect. The error probably came from 
the numerical calculation of the subprocess cross section since the the square of the 
matrix element is highly singular at such a high energy. It would not affect the cross 
section from hadron collision at Tevatron or LHC energies because the contributions 
from giuons with such a high energy center-of-mass are negligible. Therefore, we 
are confidijnt that the results of the origirial ord er-o'^ PQCD calculation by Chang 
and ChenH and the more recent oncaiiaiJ'iyiij are now in agreement. 

An alternative way to calculate the hadronic production is to use the fragmen- 
tation approximation. From general factorization theorems it is clear that, for 
sufficiently large transverse momentum Pt of the Be or i?*, hadronic production is 
dominated by fragmentation. The calculation can the n h e considerably simplified 
using this approximation, as was first done by Cheun^l3. Subsequently, the com- 
parison between the full caJc iila tio p a nd the fragmentation approximation has 



been discussed by several authorsll^'tSllj. However, very different conclusions have 
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been drawn, although their numerical results are similar. Chang et a/Jiij found that 
when \/I and px are small the difference between_the fragmentation approximation 
and the calculation is large. Berezhony et aZJU claimed that the fragmentation 
approximation breaks down even for very large Pt by examining the ratio of B* 
to Be production. Kolodziej et oLIIB claimed that the fragmentation approximation 
works well if Pt exceeds about 5 — 10 GeV, which is comparable to the Be mass, 
by investigating the Pt distribution of the Be meson. 

It is nontrivial to clarify these points since the full calculation is quite com- 
plicated, with the dominant subprocess involving 36 Feynman diagrams, but the 
importance of investigating this issue is twofold. From the theoretical perspective, 
it provides an ideal example to quantitatively examine how well the fragmenta- 
tion approximation works for calculating the hadronic production of heavy flavored 
mesons. In this process both the full a* contributions and the fragmentation ap- 
proximation can each be calculated reliably. Experimentally, it is important to have 
a better understanding of the production of the Be mesons in the small Pt region 
where the Be production cross section is the largest. The Pt distribution decreases 
very rapidly as Pt increases. 

Recently, Chang, Chen and Oake^ carried out a detailed comparative study 
of the fragmentation approximation and the full QCD perturbative calculation. 
We first studied the subprocesses carefully. By analyzing the singularities appearing 
in the amplitudes and in the Pt distributions of the subprocesses, we gained some 
important insight into the processes. We then calculated the entire hadronic pro- 
duction cross section for pp (or pp) collisions at the Tevatron energy. To facilitate 
a quantitative comparison we found it very instructive to examine the distribution 
in variable z, defined to be twice the energy fraction carried by the Be meson in the 
center of mass of the subprocess, which is in principle observable. Investigating the 
z distribution we found that the fragmentation contribution dominates when and 
only when Pt ^ Mb^ ■ It is insufhcient to evaluate the validity of the fragmentation 
approximation by only examining the Pt distribution. 

In order to obtain some insight into the process, we first focus the discussion on 
the subprocess. At the lowest order, a^, there are 36 Feynman diagrams responsible 
for the dominant gluon fusion subprocess g + g^Be + b + c. When the energy in 
the center of mass system, VI, is much larger than the heavy quark mass, the main 
contributions to the cross section come from the kinematic region where certain of 
the amplitudes in the matrix element are nearly singular; i.e., some of the quark lines 
or gluon lines are nearly on-shell. In the large Pt region, the only singularity is the 
collinear one for which the Be is collinear with fragmenting b antiquark or coUinear 
with the fragmenting c quark. Therefore, at sufficiently large Pt, the subprocess is 
dominated by fragmentation. However, when the Pt of the Be meson is small the 
produced Be, as well as the b and the c quarks, can be soft or collinear with the beam. 
In this region the square of the matrix element is highly singular because two or more 
of the internal quarks or gluons in certain Feynman diagrams can simultaneously be 
nearly on-mass-shell. Although this region is a smaller part of the phase space, these 
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nearly singular Feynman diagrams, in fact, make a large contribution to the cross 
section and dominate the small Pt region. The contributions to the cross sections 
can therefore be decomposed into two components: the fragmentation contribution, 
which dominates in the large Pt region and the non- fragmentation component, 
which dominates at smaller Pt- The contributions of these two components are 
quite clearly distinguishable in the Pt distribution of the subprocess, particularly at 
large VS. In Fig. 1 we show the Pt distribution of the subprocess when -s/I = 200 
GeV for both the full calculation and the fragmentation approximation with 
as = 0.2, TUb = 4.9 GeV, rUc = 1.5 GeV, and fs, = 480 MeV. It is easy to see 
in Fig. 1 that when Pt is larger than about 30 GeV for the Be and about 40 
GeV for the B*, the fragmentation approximation is close to the full calculation. 
However, when the Pt is smaller than these values, the deviation between the 
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Fig. 1. The Pt distributions of the Be and B* meson for the subprocess with \/l = 200 GeV. The 
solid and the doted lines correspond to the full calculation and the fragmentation approxima- 
tion, respectively. 



fragmentation approximation and the full calculation becomes large and the non- 
fragmentation component clearly dominates the production. This critical value of 
Pt is certainly much larger than the heavy quark masses, or the Be meson mass; 
it slowly increases with increasing s, which may indicate that there is an additional 
enhancement due to logarithmic terms such as Ins/m^ in the non- fragmentation 
component compared to the fragmentation component. When \/§ is not very large 
this two comp onent de composition is less distinct, since the higher twist terms can 
not be ignoredllj'EjMjlnihis case, the fragmentation approximation is not a very 
good approximationt3IU'll3. 

Let us now turn to the calculation of Be and B* production at pp colliders, 
particularly the Fermilab Tevatron. In Fig. 2, the Pt distributions of the B^ and the 
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B* mesons coming from the full calculation is compared with the fragmentation 
approxinialion. In these calculations wc used the CTEQ3M parton distribution 
functional^. For simplicity, we use a uniform choice of energy scales, choosing the 
same scale as set in the b fragmentation component; i.e. aKPr) x al{2mc). From 
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Fig. 2. The Pt distributions for Be and B* meson production at the Tevatron energy y/s = 1.8 
TeV. The solid and the doted lines correspond to the full calculation and the fragmentation 
approximation, respectively. 

Fig. 2, we see that for the meson the Pt distributions for Pt > 5 GeV are 
very similar for the full calculations and for the fragmentation approximation. 
This critical value of Pt is much smaller than that found above in the study of 
the subprocess. However, for the B* meson, the result predicted by the full 
calculation in Fig. 2 is a factor 1.5 to 2.0 greater than the fragmentation calculation 
over a much larger range of Pt, more consistent with what was found in the study 
of the subprocess. Therefore, it is difhcult to decide from only the Pt distributions 
where the fragmentation approximation is reliable for the hadronic production of 
the Be and B* mesons. The use of Pt distribution alone can be misleading. 
To clarify this issue, we introduce the distribution 

C{z)^^^ I dxidx2g{xi,fi)g{x2,n) ^ (6) 

dz J dz 

where z = 2(fci + /C2) • p/s and g{xi,fj,) is the gluon distribution function. The 
distribution C{z) provides a sensitive means to investigate the dynamics of the pro- 
duction process and the fragmentation approximation. Clearly, if the fragmentation 
approximation is valid, da{Vl, ^)/dz can be factorized as 

da(\/§,ii) sr^ ^ ^ , ^ /„x 



8 Production of the Be Meson 



where Dq.^b^{z, ^i) are the usual fragmentation functions and o-gg^g.Q. is the 
gluon fusion subprocess cross section for production of the heavy quark pair QiQi. 
In this approximation, the integrals over x\ and X2 can be performed, the frag- 
mentation function can be factored out, and C{z) is simply proportional to a sum 
of the usual fragmentation functions which is insensitive to the parton distribution 
functions and to the kinematic cuts. However, if the distribution C{z) is quite dif- 
ferent from the fragmentation functions, it is an indication that the fragmentation 
approximation is not valid. Therefore, comparing C{z) calculated in the fragmen- 
tation approximation with the full order calculation, provides a quantitative 
criterion to judge the validity of the fragmentation approximation. 

In Fig 3, the distribution C{z) calculated in the fragmentation approximation is 
compared with the full order calculation for Be and B* meson production with 
a cut of Pt > 10 GeV (Fig. 3a) and also with cuts of Pt > 20 GeV and Pt > 30 
GeV (Fig. 3b). The z distribution C{z) is sensitive to the smallest Pt region for a 




Fig. 3. The 2 distributions C{z) of the Be and B* at the Tcvatron energy y/s = 1.8 TeV. The 
solid Hnes are the full a* calculation and the doted lines are the fragmentation approximation (a) 
with the cut Pt > 10 GeV and (b) with the cuts Pt > 20 GeV and Pt > 30 GeV. 

given Pt cut because the Pt distributions of the Be and B* mesons decrease very 

rapidly with increasing Pt- From Fig. 3, some general features are evident. For 
the Be meson, with smaller Pt cuts, the distribution C{z) is overestimated in the 
higher z region by the fragmentation approximation, while it is underestimated in 
the lower z region, however, after integration over z. the result is similar to the 
full a* calculation. For the B* meson, even for the largest Pt cut, the distribution 
C (z) is underestimated the fragmentation approximation at all values of z and, after 
integration over z, the result is definitely smaller than the full a* calculation. This 
feature explains why the Pt distributions shown in Fig. 2 are similar for the Be 
meson even down to Pt ~ Mq^ but are different for the B* meson. This also shows 



Production of the Be Meson 9 



that it is simply fortuitous that the Pt distribution of the Be calculated in the 
fragmentation approximation is similar to that from the full calculation for Pt 
below this value, particularly down to Pt ~ Mb^- It is also clear that when Pt is 
increased the distributions become closer. As shown in Fig. 3b, when Pt is as large 
as 30 GeV the curves calculated in the fragmentation approximation are quite close 
to the full calculation. This indicates that the fragmentation approximation is 
valid in the large Pt region, as expected. We emphasize here that the difference 
between the fragmentation approximation and full calculation is not universal, but 
is process-dependent. 

Finally, we examine the ratio of cross sections for B* production calculated in 
the fragmentation approximation with the results of the full af calculation. As 
discussed above, for B* meson production the fragmentation approximation always 
underestimates the full result. The deviation from the full calculation for the B* 
meson can be used as a criterion to test the validity of the fragmentation approxi- 
mation. The results for the total cross section (t{Pt > PTmin) for various Pt cuts 
are listed in Table I. Taking agreement within 30% as the criterion for the validity 
of the fragmentation approximation we also see that Pt should exceed about 30 
GeV, a value considerably larger than the heavy quark masses. We note that this 
conclusion is rather insensitive to the choice of the energy scale fi and the parton 
distribution functions. 

Table 1. Total cross sections u{Pj'b^ > PTmin) in for hadronic production of the Be and 
the B* mesons predicted by the calculation and the fragmentation approximation assuming 
various Pt cuts and \Y\ < 1.5. The CTEQ3M parton distribution functions were used and the 
values /s^ = 480 MeV, = 1.5 GeV, mt = 4.9 GeV, and Mb^ = 6.4 GeV. 



Prm»n (GeV) 



"B^ifrag. ) 
o"i3; ( frag. ) 

g-g; ( frag. ) 

^B'Jai) 







10 



15 



20 



30 



1.8 0.57 0.087 0.018 4.8 x 10" 



1.4 
4.4 
2.3 



0.47 0.071 0.014 4.0 x 10" 



1.4 
0.78 



0.22 
0.12 



0.52 0.55 0.56 



0.041 
0.025 

0.61 



6.3 X 10"* 
5.3 X 10-* 
1.1 X 10-2 1.3 X 10-3 
6.8 X 10-3 9.2 X 10-* 



0.63 



0.70 



In summary, from the study of both the parton subprocess and the hadronic 
process, we can conclude that the fragmentation mechanism dominates when and 
only when Pt ^ Mb^ ■ This conclusion is independent of the choice of the energy 
scales and the parton distribution functions. It is only fortuitous that the Pt 
distribution of the Be in the fragmentation approximation is similar to that of the 
full calculation for Pt as low as Mb^- 



4. Conclusion 

Both the spec troscopy o f be system and the decays of the Be mesons have been 
widely studiedllB£3'Eil'E3. The excited states below the threshold will decay to the 
ground state Be by emitting the photon(s) or n mesons. The golden channel to 
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detect the Be meson is Be — » J/ip + i^{p)- However, the branching ratio is quite 
small, around 0.2. The exclusive semilcptonic decay mode Be — *■ J /ip -\- I + vi 
has a relatively larger branch ratio, but there is "missing energy". Comparative 
studies indicate that while fragmentation approximation works very well for the 
production of the Be and B* in high energy e^e~ collisions, the condition for the 
applicability of the fragmentation approximation in hadronic collisions is that the 
transverse momentum Pt of the produced Be and B* be much larger than the mass 
of the Be meson; i.e., Pt ^ Mb^- Theoretical studies for the production of the 
Be and B* mesons indicate that an observable number of Be and B* events can be 
produced at LEP I and at the Tevatron. With the discovery of the Be in the near 
future, a new page in heavy flavor physics will be opened and the study of the Be 
meson will no longer be academic. 
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